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We study heat conduction behavior of one-dimensional lattices with asymmetric, momentum con- 
serving interparticle interactions. We find that with a certain degree of interaction asymmetry, the 
heat conductivity measured in nonequilibrium stationary states converges in the thermodynamical 
limit. Our analysis suggests that the mass gradient resulting from asymmetric interactions may 
provide a phonon scattering mechanism in addition to that caused by nonlinear interactions. 
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The heat transport properties of low-dimensional sys- 
tems have attracted intensive studies for decades |ll420l| 
(see also Refs. j2l| - |23| for reviews and references therein). 
A challenge is to relate the heat conduction behavior of 
a system to its microscopic ingredients. In 1984 Casati 
et al. investigated the role chaos may play [4] , and since 
then their seminal work has trigged numerous efforts for 
identifying the microscopic mechanism(s) of the Fourier 
law. In a one-dimensional (ID) case, the Fourier law 
states 



J 



dT 
dx ' 



(1) 



where J is the heat current, ^ is the spatial tem- 
perature gradient, and k is a finite constant termed as 
"thermal conductivity." The heat conduction behavior is 
also known as " normal heat conduction" if it follows the 
Fourier law or "abnormal heat conduction" otherwise. 
Now it has been clarified that chaos is neither sufficient 
nor necessary to the Fourier law [5|-|7|. 

For ID lattices, another significant step was made in 
1998 by Hu et al., who pointed out that, besides the 
dynamical properties, whether or not the system has 
a conserved total momentum is another key ingredient 
[9|, |lO|; i.e., lattices with (without) a momentum conser- 
vation property should disobey (obey) the Fourier law. 
In 2000 Prosen and Campbell went a step further; they 
proved that for ID momentum conserving lattices with 
non- vanishing internal pressure the heat conductivity di- 
verges in the thermodynamical limit |8J. Though for 
lattices with a vanishing internal pressure their proof 
is not applicable, many numerical studies support the 
same conclusion. In addition, in their later study Prosen 
and Campbell also showed that momentum conserving 
is not a necessary condition for abnormal heat conduc- 
tion [7|. More recent progress was made by employing 
the renormalization group analysis for hydrodynamical 
models fly, [i3| and the mode coupling theory [18l - l2l| . 
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Again both theories predict a divergent heat conductiv- 
ity in ID momentum conserving systems. This progress 
has greatly deepened our understanding of the heat con- 
duction problem. However, in spite of this fact, to our 
knowledge there are also three counterexamples jlll - [l3l | 
which are momentum conserving, but instead have con- 
vergent heat conductivity. (Reference [1^ , closely related 
to the present study, will be discussed in detail later.) 

In this paper we investigate the effects of asymmetric 
interparticle interactions on the heat conduction of ID 
momentum conserving lattices. It has been well known 
that asymmetric interactions are important for lattice 
systems; e.g., they may induce a non- vanishing inter- 
nal pressure and thus the thermal expansion effect [2j|. 
Asymmetric interactions may have significant implica- 
tions on transport properties as well. This was implied 
by the studies in Ref. [8| and later was well shown in 
Refs.|18l.[l9j. where two generality classes, corresponding 
to whether or not the interactions are symmetric, were 
put forth. In Ref. [1^1, by proposing a "mode cascad- 
ing" relation and incorporating it into the mode coupling 
theory, the authors theoretically predicted and numeri- 
cally verified that for systems of symmetric interactions, 
or equivalently, for systems with equal specific heat ca- 
pacities at fixed pressure and volume, the bulk viscosity 
is finite in the thermodynamic limit while for those of 
asymmetric interactions it is divergent. In addition, the 
heat conductivity is suggested to diverge in both cases 
but in different ways. Another development of the mode 
coupling theory, i.e., the self-consistent mode coupling 
theory [l^, has led to a the consistent conclusion. In 
particular, it concluded that for the ID momentum con- 
serving systems the heat conductivity would diverge and 
the divergent exponent is i or i respectively, for sys- 
tems with leading cubic or quartic anharmonic nearest 
neighbor potentials. 

In contrast to these results, in the following we shall 
show that in general momentum conservation does not 
necessarily imply the breakdown of the Fourier law. Our 
key finding is the existence of the converged, finite heat 
conductivity in ID lattices with asymmetric interparti- 
cle interactions. (In the following we refer to "lattices 
with asymmetric interparticle interactions" as "LWAH" 
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FIG. 1: (Color online) The schematic plot of the potential 
function V{x) given by Eq. ((3]) for r > 0, r — 0, and r < 0, 
respectively. 
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for short.) We shall present our simulation results first, 
then discuss their relation to existing theoretical and nu- 
merical studies. 

We consider homogeneous lattices with nearest neigh- 
boring coupling, whose Hamiltonian is 



FIG. 2: (Color online) The heat conductivity k vs the num- 
ber of particles TV in our lattice model for various values of 
the interaction asymmetry parameter r. The size and tem- 
peratures of the two heat baths coupled to the system are 
Na = 12, Tt = 3, and Tr = 2, respectively. The error bars 
(not shown) are much smaller than the symbols. The dashed 
line indicates ~ A''. 
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where p^ and xi are, respectively, the momentum and po- 
sition of the ith particle, and V is the potential for the 
interparticle interaction. As no on-site potentials are in- 
volved, this is a momentum conserving model. We as- 
sume that the component particles are identical and have 
a unit mass /i, and the lattice constant, denoted by a, is 
a unit as well. For our aim here the interaction poten- 
tial with an adjustable asymmetry is favorable. We have 
investigated several different forms of the interaction po- 
tential, some of which will be discussed later, but with 
all of them qualitatively the same results have been ob- 
tained. So as a typical example we shall focus on the 
following potential; 



V{x) = -(a; + r)2+e" 



(3) 



Here r is a controlling parameter that governs the de- 
gree of the interaction asymmetry; by increasing \r\ from 
zero where the potential is harmonic and symmetric, one 
gets increasingly stronger asymmetry. Fixing the system 
size to be that at zero temperature with a free boundary 
condition, the potential asymmetry implies a nonzero in- 
ternal pressure at a finite temperature: While for r > 
the internal pressure is positive and the system is ther- 
mally expansive, for r < it is negative and the system 
is of negative thermal expansion. Note that a; = is the 
equilibrium point of the potential, and Vix) for r and — r 
is symmetric with respect to a; = 0. The schematic plots 
of the potential function are presented in Fig. [1] 

To measure the heat conductivity of our system, two 
Nose-Hoover heat baths [2^ at temperatures Tl and Tr 
are coupled to the left- and rightmost A^o particles, whose 
motions follow Xi = —.pi = —jM: — 'i±Pi, and <^± = 



Pi 



j-njf 1. The Boltzmamn constant ks is set to be unity. 

The motions of TV particles between the heat baths are 



governed instead by i^ = — and pi 



- ^ . Given these 



motion equations, the evolution of the system can be 
simulated by standard numerical integrating algorithms. 
In our simulations, initially all the particles are as- 
signed to reside on their equilibrium positions with a 
given random velocity generated from the Maxwellian 
distribution at an average temperature T = \{Tl + Tn). 
Then the system is evolved for a long enough time (> 10® 
for all the cases investigated) to ensure that it has relaxed 
to the stationary state. After that the next evolution of 
time ^ 10^ is performed to obtain the time average of the 

following quantities: (i) local temperatures T^ = -^r-^; (ii) 
local heat currents Ji = (ii^j-) as adopted convention- 
ally p, [I4I ; and (iii) heat conductivity k based on 



JNa 
AT 



(4) 



by assuming the Fourier law [see Eq. ([T|)]. Here J = (Ji) 
and AT = Tl — Tr. Before we proceed, we emphasize 
that the numerical results to be presented do not depend 
on the simulation details considered here. In particular, 
we have verified that within the error range they do not 
change if the relaxing time and the average time are in- 
creased (by five times) , or if we use different definitions of 
the local heat current, e.g., Eqs. (17) and (23) in Ref. [21| 
are taken. This is also the case if the leapfrog integrating 
algorithm mainly adopted in this study is replaced by the 
Runge-Kutta algorithm of seventh to eighth order. 

Our main results are summarized in Fig. [21 where 
the dependence of k on the system size N is studied for 
various values of parameter r. The most striking fact 
revealed there is that for |r| > 1 the heat conductiv- 
ity becomes size independent for TV > 10**, suggesting 
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FIG. 3: (Color online) The temperature profiles for r = 1.5. 
The size and temperatures of the two heat baths coupled to 
the system are A'o = 12, Tl = 3, and Tr = 2, respectively. 



that the Fourier law holds. This is o ppo site to the the- 
oretical % H [H-Illl and simulation [ll, [H, ^ results 
that in ID momentum conserving lattices the Fourier law 
does not hold. To give further support for the converging 
heat conductivity observed for \r\ > 1, we plot in Fig. [3] 
the temperature profiles for r — 1.5. It shows that for 
iV > lO"* the temperature profiles can be well rescaled 
by (^). This fact justifies the calculation of the thermal 
conductivity based on Eq. (|3]) when the system size is suf- 
ficient large (T), IH, |l3] • In addition, we find that the heat 
conductivity is the same for r and — r, suggesting that in 
our model thermal expansion and negative thermal ex- 
pansion have the same implication for heat conduction. 
As a comparison the heat conductivity for the harmonic 
chain (with r = 0) is presented in Fig. [5] as well; it di- 
verges linearly with the system size, as expected [l|. We 
have also studied other asymmetric potentials and found 
qualitatively the same results. For example, in the case 
of V{x) = (1 + X)x^ for a; < and (1 - X)x^ for x > 0, 
where < |A| < 1 serves as the asymmetry controlling 
parameter, k has been observed to saturate in the inves- 
tigated parameter range 0.5 < A < 0.8 when the system 
size is large enough {N > lO'*). For this reason we con- 
jecture the finite conductivity is a general existence in 
ID LWAII. 

Given this one may wonder why previously this was 
not observed in numerical studies or predicted by theo- 
retical approaches. For the former our analysis suggests 
the reason could be that the sizes of the systems investi- 
gated in previous simulations are not large enough. The 
transient system size, denoted by N*, where for N > N* 
the heat conductivity becomes saturated, is found to de- 
pend sensitively on the asymmetry parameters. In the 
example presented in Fig. [51 we notice that N* takes 
its minimal value at \r\ « 1 and increases quickly as \r\ 
becomes smaller. Hence for a less asymmetric potential, 
e.g., r = 0.5 (see Fig. [5]), a much larger N* (> 10^) is 
expected. Moreover, for a less asymmetric potential, the 



FIG. 4: (Color online) The mass density function, respec- 
tively, for our model with r — 1.5 (dashed line) and r — —1.5 
(dotted line), FPU-/? model (solid line), and the variant ding- 
a-ling model [131 (dashed-dotted line). The system size is 
A'' = 1200. Other parameters are Tl ^ 3, Tr = 2, and 
A'^o = 1- For the variant ding-a-ling model the natural fre- 
quency of the springs connecting the even numbered particles 

is U!o — 1. 



heat conductivity seems to depend on the system size in 
a power law for N <IC N*, which explains why in pre- 
vious studies a power law divergent, rather than a con- 
vergent heat conductivity, has been found in ID LWAII. 
Indeed, as an example it is easy to check that the asym- 
metry of the Fermi-Pasta-Ulam(FPU)-a-/3 model with 
Vai3{x) = ix^ -I- ^x^ + jx^ as was considered in Ref. [l^ 
is much weaker than the case of |r| = 1 in our model. 

On the other hand, the existing theoretical predictions 
may not be applicable to the LWAII. It should be noticed 
that in these theoretical treatments, the system is usu- 
ally assumed to be at an equilibrium state with a uniform 
temperature, and thus a homogeneous mass distribution. 
But, however, in the LWAII there is an important dif- 
ference between nonequilibrium stationary states (with 
a temperature gradient) and equilibrium states: In the 
former the thermal expansion effect may simultaneously 
give rise to a mass gradient across the system. This is 
essentially different from lattices with symmetric interac- 
tions where a mass gradient is not expected in either the 
equilibrium or the nonequilibrium cases. In Fig. [31 the 
mass density function p for our model is compared with 
that of the FPU-/3 model with Vpix) = ^x^ + \x^. It 
shows clearly that, when being coupled to two heat baths 
at different temperatures, a mass gradient is eventually 
established in our system for r ^ when the stationary 
state is approached. It has been known that in systems 
with symmetric interactions, a nonlinearity of interac- 
tions may result in scattering to the heat current that 
is strong enough to establish the temperature gradient 
but not strong enough to lead to normal heat conduction 
[a, [3|- Therefore in systems with asymmetric interac- 
tions, the resultant mass gradient may provide an addi- 
tional scattering mechanism to the heat current. We con- 
jecture that this is the reason why normal heat conduc- 



tion can then be observed. According to the fluctuation- 
dissipation theorem, such a macroscopic, nonequiUbrium 
effect must have its microscopic, equiUbrium counterpart, 
but the latter may have not been the object of the exist- 
ing theoretical studies. 

As one more evidence for our conjccturc-that the mass 
gradient may provide an additional scattering mechanism 
to the heat current-it is worthwhile to notice that in a 
very recent study [lj|, a different ID momentum conserv- 
ing system having a finite thermal conductivity has also 
been reported. The system is a momentum conserving 
variant of the "ding-a-ling" model [2|: The even num- 
bered particles are bound to the adjacent even numbered 
particles by harmonic springs, and are subject to the 
elastic collisions with their neighboring odd numbered 
particles. The odd numbered particles are free except 
for elastic collisions with their even numbered neighbors. 
Significantly, the interactions are asymmetric due to the 
elastic collisions, and as a result the system is thermally 
expansive. In Fig. 4 the mass density of the system in a 
nonequilibrium state is compared with our system; it can 
be seen that its asymmetry degree of the interactions is 
even stronger than our system with \r\ = 1.5. This ex- 
plains why the saturating regime of the heat conductivity 
can be numerically accessed in this system [l3| , provided 
that our conjecture is correct. 

Finally, we would like to emphasize that, not only in 
ID LWAII but also in ID gases, the mass gradient in 
nonequilibrium stationary state may have significant ef- 
fects on heat conduction. This has been shown in a ID 
hard-core gas with alternative molecule masses [2g| and 
a variant ID hard-core gas model where two neighboring 
molecules are bound by a massless string [27]. The inter- 
particle interactions in both of them are asymmetric. In 
the former it has been shown both analytically and nu- 
merically that, when the system is exposed to two heat 
baths of different temperatures, the temperature gradi- 
ent across the system is maintained by the mass density 



gradient. In the latter the heat conduction behavior has 
been found to dramatically depend on whether or not the 
system has a nonzero external pressure (equivalently, a 
non-zero internal pressure due to the force balance). In 
addition, as was stressed in Ref. j27| , the heat transport 
properties measured in equilibrium and nonequilibrium 
states could be qualitatively different. However, it should 
be noticed that in these two studies the heat conductivity 
has been shown to diverge in the thermodynamical limit, 
which implies that, lacking a phonon scattering mecha- 
nism in gases, the mass gradient cannot guarantee the 
Fourier law exclusively. 

To summarize, we have performed a numerical investi- 
gation for several ID momentum conserving LWAII and 
observed normal heat conduction behavior. Comparing 
our finding in ID LWAII and the heat conduction charac- 
teristics of ID lattices with symmetric interactions such 
as the FPU-/? model, we conjecture that the mass gra- 
dient may provide a phonon scattering mechanism in 
addition to that which is caused by nonlinear interac- 
tions, which jointly leads to the observed normal heat 
conduction behavior in ID LWAII. Based on our under- 
standing, we conjecture the same mechanism also works 
in the two-dimensional (2D) case. Indeed, normal heat 
conduction has been observed in both our ongoing study 
of 2D momentum conserving LWAII [28] and a recent 
numerical study of thermal conductivity in empty and 
water-filled carbon nanotubes [2^. As thermal expan- 
sion is ubiquitous among real lattice systems, suggesting 
their interparticle interactions are generally asymmetric, 
we expect that the Fourier law generally holds in real 
low-dimensional systems. In this regard experimental in- 
vestigations of carbon nanotubes and graphene flakes of 
large sizes are very desirable. 
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